Exact ground states for a class of one-dimensional frustrated quantum spin models 
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We have found the exact ground state for two frustrated quantum spin- 1/2 models on a linear 
chain. The first model describes ferromagnet - antiferromagnet transition point. The singlet state 
at this point has double-spiral ordering. The second model is equivalent to special case of the spin- 
1/2 ladder. It has non-degenerate singlet ground state with exponentially decaying spin correlations 
and there is an energy gap. The exact ground state wave function of these models is presented in a 
special recurrent form and recurrence technics of expectation value calculations is developed. 
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I. INTRODUCTION 



Last decade frustrated Heisenberg models have been a subject of intensive studies 1-12 . Of main interests are ground 
state properties with respect to variations of exchange integrals and character of the phase transitions. In particular, 
these properties may be important in the theory of high-T c superconductivity. 
^ \ There is much interest in quantum spin systems with competing interactions for which exact ground state can be 
constructed. The first example of such a model has been given by Majumdar and Ghosh 13 . They have considered 
s = h chain with antiferromagnetic nearest- and next-nearest neighbor interactions and the strength of the second 
interaction is one half of the first. The ground state of this model is two-fold degenerate, consists of dimerized singlets 
and there is a gap in the spectrum of excited states. 

Another example found by Affleck, Kennedy, Lieb and Tasaki is an S = 1 spin chain with special bilinear and 
biquadratic interactions (AKLT model). This model has unique resonating- valence-bond ground state and ground 
state correlations have exponential decay. Besides, there is the gap between ground state and excited states. Further 
Q ■ generalizations of the AKLT model have been studied in a number recent papers 15 . 

In this paper we present two classes of s — i chains for which exact ground state wave function has a special 
recurrent form. These models have competing ferro- and antiferromagnetic interactions and their ground states can 
' be either ferromagnetic or singlet depending on the relation between exchange integrals. The first type of exactly 
. solvable models is related to the systems at the ferromagnet-antiferromagnet transition point when the ferromagnetic 
and the singlet states are degenerate. The calculation of the spin correlation function in the singlet ground state 
^ — I \ shows the spiral magnetic order at this point. 

The model of the second type has the nearest, next-nearest and next-next-nearest neighbor interactions depending 
on one parameter. This model has the non-degenerate singlet ground state for cyclic chains and for a certain region 
of the parameter and its ground state properties are similar to that of AKLT model. In other words, this model has 
' all properties of the "Haldane scenario" 16 , though it is the model with half-integer spin. We note, however, that the 
considered model has two sites in an unit cell and it is equivalent to the special case of a spin ladder. In some limit 
this model reduces to the effective spin-1 chain for which the ground state wave function coincides with that for the 
AKLT model. 

The paper is organized as follows. In Section 2 we will consider the model of the first type and describe the exact 
singlet ground state wave function as well as details of the spin correlation function calculations. Section 3 and 
Appendix present the study of the model of the second type. The results of the paper will be summarized in Section 
4. 



II. FRUSTRATED SPIN CHAIN AT FERROMAGNET-ANTIFERROMAGNET TRANSITION POINTS. 



A. The exact ground state wave function. 



Let us consider s — i spin model with nearest- and next-nearest neighbor interactions given by the Hamiltonian 

M M ^ N 

H=- X)(S ai _i • S 24 - -) + J 23 £)(S ai • S 2m - -) + J 13 £(S, • S 4+2 - -), (1) 

i—1 i—1 i—1 

with periodic boundary conditions and even N — 2AI. 
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FIG. 1. T = phase diagram of the model (1). The solid line is the boundary between the ferromagnetic and singlet phases. 
Model (1) along this line is given by Hamiltonian (2). Circles correspond to the special cases of this Hamiltonian. 



If J 23 < 1) then the ground state of (1) is ferromagnetic (singlet) at 5 < (5 > 0), where S = Ji 3 + 2 (i- 2 ,/ 23 ) (Fig-1)- 
The equation 6 = defines the line of transition points from the ferromagnetic to the singlet state, when energies 
of these states are zero. The model (1) along this line is given by the Hamiltonian depending on the parameter v 
{v > 0): 

M M _ N 

H=- 5>2i-i ' S 2 « ~-)-{y-\) ]T(S 2l • S 2l+1 --) + — J2(s l ■ S l+2 -), (2) 

i—1 i—1 i—1 

with periodic boundary conditions. 

We note that the Hamiltonian (2) has a symmetry: its spectrum coincides with the spectrum of H(y) obtained by 
the following transformation 

H{v) = {y - 1) ff(^y), ">1 

This transformation permutes the factors at the first and the second terms in the Hamiltonian (2). Thus, due to 
the symmetry it is sufficient to consider the range < v < 2. 

First, we will show that the ground state energy of (2) is zero. Let us represent the Hamiltonian (2) as a sum of 
Hamiltonians H n of cells containing three sites 

M 

H = J2(H2i-i+H 2i ), (3) 

i=l 



where 



H 2 i-i - --(S 2 j-i • S 2j - -) — (S 2j • S 2j+ i - -) + 2 ^ (S 2 »-i • S 2j+ i - -) 

H 2i = -^(Sli+l ' S 2i+2 - — (S 2 « ' S 2i+ i - -) + 2 ^ (S 2 » • S 2 j+2 - 



Eigenvalues of each H n are 



v 2 - v + 1 

Ai = A 2 = 0, A 3 = > 

2v 
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We will present a singlet wave function which is the exact one of each H n with zero energy and, therefore, it is the 
exact ground state wave function of (2). This function has a form 

*o(M) = P *m, (4) 

where 

*Af = {sf + Vs\ + ZASg . . . + ^s^)(s^ + vsj .. . + issjj) . . . 

x(4„-i + ^4„--- + ^4)---(4-i + ^4) I 1.2....JV), (5) 

where sf is the s = \ raising operator. 

Eq.(5) contains M = ^ operator multipliers and the vacuum state | 1,2,... N) is the state with all spins pointing 
down. The function \I>m is the eigenfunction of S z with S z = but it is not the eigenfunction of S 2 . P is a projector 
onto the singlet state. This operator is 17 



p o = t, — o / da / d(3 / sin 7(^7 
8?r 2 Jo Jo Jo 



e iaS z e i 1 s Xe i0S z 



where S x (z) are components of the total spin operator. 

The function ^ m contains components with all possible values of spin S (0 < S < M) and, in fact, a fraction of 
the singlet is exponentially small at large N. This component is filtered out by the operator P . 

It is not difficult to check that 



H n ^ M = (6) 

for n = 1, ... (N — 2) and, therefore, the ground state energies for all values of spin S of an open chain describing by 
the Hamiltonian 

N-2 

H op = ^ H n (7) 

n=l 

are zero. 

The operators i?jv-i and Hn do not give zero acting on ^ m but 

H N _ 1 ( N - ) P ' ( 5>M = -Po^AT-l(AT)^'Af = (8) 

The latter equation can be easily checked using the fact the first bracket in (5) can be replaced by (1 — v)sf under 
the projector Po- 

Eqs.(6) and (8) mean that \&o(M) is the exact singlet ground state wave function of (2) for any M and the ground 
state energy is zero. We note that the ground state energy coincides with its exact lower bound because H is the sum 
of non- negative defined operators (Eq.(3)). Of course, the trivial ferromagnetic state has zero energy as well. 

In particular case, v = 2, when J\i = J 2 3 = —1 and J13 = j, another form of the exact singlet ground state wave 
function has been found in Ref.18. It reads 

* = 2M[M[m,n]... > 

where denotes the singlet pair and the summation is made for any combination of spin sites under the condition 
that i < j,k < l,m < n . . . . However, it is not clear how this function can be generalized to v 7^ 2. 
The following general statements relevant to the Hamiltonian (2) can be proved: 

1) . The ground states of open chains described by (7) in the sector with fixed total spin S are non-degenerate and 
their energies are zero. 

2) . For cyclic chains the ground state in the 5* = sector is non-degenerate and has momentum it. The ground 
state energies for < S < M are non-zero. 

3) . The singlet ground state wave function for open and cyclic chains coincide with each other. 

4) . The singlet ground state wave function ^o(M) is superstable 19 with respect to any cell operator H n , i.e. ^/o{M) 
is the ground state wave function of the Hamiltonian H + \H n for —1 < A < 00. 
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B. The norm of the ground state wave function. 



Let us return to the problem of the projection of the function \I>m- As one can see from Eq.(5) the function 
satisfies a recurrent equation 

*m = [(4 + "4) + » S+ ( N ~ 2 )] I 1. 2)*m-i, (9) 

where ^m-i is the function (5) for the system of (N — 2) spins on sites 3, 4, . . . N and S + (N — I) = J2iLi+i s t ■ 
In principle, it is possible to generate ^q(M) starting from M = 1 and using Eqs.(4) and (9). However, it is more 

convenient to obtain the recurrent formulae for expectation values (norm and correlators) with respect to the function 

^o(M) directly. 

First, we consider a norm of ty (M) which has a form 

G M = (*o(M)* (M)) = (* A fP *M> = \ I $ M (7)sin 7 d7, (10) 

* Jo 

where 

$M(7) = $M(7l>72) l 7 l=72=7 

and 

$m(7i,72) - (*mcxp^S + (A0+z^S-(A0)*m> ( n ) 
Commuting operators in Eq.(ll) and using the fact, that S z ^m — 0, we rewrite Eq.(ll) in a form 

$m(7i,72) = (y M cM^S-(N)) C xp(iz'S+(N))<f M ), (12) 

where 



F/2. V7172 , _ pfism y/Yi 72 

Z — a j tcHl — . Z — 



7i 1 V 72 ^ 

It follows from Eqs.(9) and (12) that the function <&m(7i, 7 2) satisfies the equation 

$m(7i,72) = [1 + v 2 (1 + vfzz' + v(y + 1)(1 - zz')(z^ + z' A) 

d 2 

-^(l-zz 1 ) 2 ^-^} $ M -i(7i,72) 
Using this equation we obtain for $m( 7 ) 

$m(7) = [1 + v 1 - (1 + vf sin 2 \ + v{v + 1) sin 7 -f- 

_^ CQS 2 + J_d _ 

2 d~f 2 sin 7^7 

The solution of Eq.(13) is 

$ M (y) = i M (y)$o(y), (14) 

where 

L(y) - ^ + - v{v + i)(i- - y (i + .) 2 ^(i - y)^, (15) 

y = cos 7 and <&o (y) = 1. 

According to Eq.(15) $m(v) is a polynomial in y of order M. It turns out that further calculations will be simplified 
if (y) is expanded over Legendre polynomials P n (y) : 

M 

$ M (y) = 5>„(M)P„(y) (16) 

n=0 
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The coefficients c n (l) are defined by the recurrent equation 

„ .,, n {vn + lf .,. ^ 2 (n 2 + n) + (y - l) 2 .,. n + I {vn + v - if .,. , . 

c„(/ + 1) = 2n _/ 2 C-i(Q + ^ ^ L cnd) + — 3 [ - 2-^^(0 (17) 

with initial condition co(0) = 1. Besides, c„(i) = at n > I. 
The norm Gm is given by 



G M = 7 : <$> M {y)dy = c (M) (18) 



C. Spin correlations. 

The spin correlation functions can be found in the same way as Gm- It is convenient to express the scalar product 
Si • Sj by the permutation operator pij = 2Si ■ Sj + 5. Then the non-normalized expectation value Q(i,j) = 
(*o(M)(S< • Sj - |)*o(M)) is 

Q(i.i) = -\(MM)(p id - l) 2 *o(M)) = -i(* M ( Pi;j - l)P (p 4J - 1)* M ) (19) 

The equations for Q(i,j) can be obtained as in the derivation of Eqs.(14)-(15). They are somewhat different for 
even and odd i. For example, the equations for Q(l, n) have forms 

0(1, 21) = -^=^- f\ l -±±{L l -- Ll L M - 1 l)dy, (20) 



Q{l,2l-l) = - v2{ \ 1)2 J^llVtf-'LtL"-* l)dy, (21) 



where 



1-2*1-,) <-,>(l +,)<<!-,,)<. 



It is clear, that 

Q(2,n + 2) = Q(l,N-n + l) 
Making use consequent integration of Eqs.(20) and (21), we obtain 



|* i±^(L'- 2 L 1(2) i M -'l)rf y = J\ l±l { L l -H){L m L M - l l)dy = f {L 1(2) L l - 2 l){L M - l l)dy, 



where 

= 2 ' 2 y " v "" i/vi y 2 y ^ ' dy y ^ "'dy 1 



L = Li = - + 1 S - ^ - 1)(1 - y 2 )-r, ~ -=-(1 + y) 2 3-(l - I/)-j-> 



7 1 + ^ n 2^ d (1 + y) 2 d , d 



So, Eqs.(20) and (21) can be rewritten as 

0(1,20 = - ( "~ 1)2 y $1-1(1/) $M-;(y)^, (22) 
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Q(l,2Z-l) = - ^ 4 1)2 J (L 2 *,_ 2 (y)) *M-j(l/)d», (23) 

where 

*,(») = (24) 

The function $/ (y) can be expanded over Legendre polynomials similarly to <I>; (y) 

i 

= $>n(0^n(i/), (25) 
n=0 

and coefhcients a n (l) satisfy the equations 

,, n (vn + v-1) 2 v 2 (n 2 + n) + 1 ,,, n + 1 (i^n+1) 2 
° n( + } = 2^1 2 -On-i(0 + — ^ MO + ^-«n+i(0 (26) 

with initial condition oo(0) = 1 and a n (l) = at n > I. 

Using Eqs.(18), (22) and (23) we can express the correlation function K(l,n) = (Si • S„)/Gm in the forms 

K(l, 21) = 1 - ^ V (27) 
V ' ' 4 2c M ^ 2n + l V 1 



2Z " 1} - I ~ 2co (M) 2. 2^T1 ' (28) 

v ' n=0 

where coefficients b n (l — 1) are defined by a n (l — 2) as follows: 

- 1) = ™ 70n-i(i - 2) + + 9 + a w (/ - 2) + i^^ an+l( ; - 2) (29) 

2 2n — 1 2 2 2n + 3 

Therefore, the calculation of the spin correlation function reduces to the solution of the recurrent equations (17) 
and (26) which were used for numerical calculations of the spin correlation function for finite systems. 
At large M the solutions of the recurrent equations (17) and (26) have scaling forms: 

c„(M) = 2M(M!)V M sexp f Mf (s) + — - InM + f^s) J , (30) 



a n (M) = 2M(M!)V M scxp (m/ (s) - - InM + gi (s)j , (31) 

where the parameter 

_ 2n+ 1 
S ~ 2M 

can be considered as a continuous variable. 

We note that Eqs.(30) and (31) are not valid for special values of v, v = (m = 0, 1, 2 . . .). For these v the last 
term in (17) vanishes when n — m and Eq.(17) reduces to (m + 1) equations for c n (l) with n < m. 

The recurrent equations (17) and (26) at M » 1 and v ^ reduce to differential ones for fo{s), fi(s) and 
<7i(s). For example, /o(s) satisfies the equation 

( t i x df (s)\ 2 2 {ldf (s)\ 
exp (/„(,) -,_j= a cosh 

with initial condition /o(l) = — In 4. 
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Its implicit solution is 



/ (s) = -21n£ + 2^1ntan|, (32) 



where 

sin£ 

As it follows from Eq.(32), a n (M) and c„(M) as a functions of n have a sharp maximum at -jg 
The functions fi(s) and gi{s) are given by 



/i(s) = — lnf + -^lnsinf - -ln(-|2 — ) + C— — , (33) 



/ \ ^ — xJis , . Of — 4: , . J. _ . Sill C tua L . f — J_ , . 

5i (*) - — InC + In sin ^ - - In(-^ - -^) + C— — , (34) 



2 — 3^, . Zv — 4 , . . 1 , ,sin£ cos£. „v—\ 
' -^ lnsm ^2 ln( — " — 
where constant C ~ 0.969. 

To obtain K(l, 21) at N — > oo we substitute Eqs.(30)-(34) into Eq.(27) and replace the sum over n by the integral 
over s. This integral is calculated by the method of steepest descent. The saddle point is so = ^sin(^) and the 
integrand does not depend on parameter v. The final result for the spin correlation function K(l,2l) at N — > oo is 
also independent on v. 

A - (1 ,^H^)+0(1) (35, 
For the particular case v — 2 Eq.(35) has been obtained earlier in Rcf.18. 

The corresponding calculations for K(l, 21 — 1) to within terms ~ lead to the same expression (35). But taking 
into account terms ~ ^ we find that the difference K(l, 21 — 1) — K(l, 21) is 

A-,l,2 ( -l)-A'( 1 ,2 i )^^ 8 in(f) +0 (i 5 

The latter equation means that the double-spiral structure exists. The pitch angle of each spiral is ^ and there is 
a small shift angle Ai^s = j^^r between them: 

Jr(l,20 = ^co S ( 27r(2 ^~ 1) -Ay), (36) 

K(l,2l + l) = \cos(^p) (37) 
This shift angle reflects the fact that the unit cell contains two sites unless v = 2. 

Eqs.(36) and (37) show that the long range spiral order exists in the singlet ground state of the Hamiltonian (2) in 
the thermodynamic limit and the double-spiral state is formed. 

It is interesting to note that correlators (36) and (37) coincide with those obtained by using the simple 'quasi- 
classical' trial wave function in a form 

V> c ; =cxp (^2^n4^J \1,2,...N), 

where 

/27m (-l) n nv- 
<T„ = exp -— + 



N N v 

Thus, the quantum ground state of the large- iV limit resembles the classical one though for small size systems the 
quantum fluctuations are essential. 

The formation of the spirals having the period which is equal to the system size reflects the tendency to the creation 
of the incommensurate spiral state at the antifcrromagnctic region when 5 > 0. The behavior of the system in the 
vicinity of the transition point has been studied by us 7 for the model (1) with J\i = J23 = — 1, J13 = j + S. For S <C 1 

the period of the spiral is finite and is proportional to . The transition from the ferromagnetic to the singlet state 
is a phase transition of the second order with respect to S. 
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D. Special cases of the model. 



There are the special points, v = ^ (m = 1, 2, . . . ) at which Eqs.(36) and (37) are not valid. 
At v = \ the function (4) reduces to the product of singlets 

tt (M) = [2,3][4,5]...[JV,l], (38) 

and K(2, 3) = K(A, 5) = . . . = K(N, 1) = — |. Other spin correlators are zero. 

Analysis of Eqs.(17), (26) and (27) shows that the ground state correlations of the model (2) with v = ^ (m > 3) 
have antiferromagnetic character with an exponentially decay: 

71 

tf(l,n+l)~(-l) n exp( ), 

where the correlation length is 

r e = -21n- 1 (l- 2 ) (39) 
\ m(m — 1) / 

The crossover between the spiral state at ^ < v < -^zi an d the antiferromagnetic state at v = ^ occurs in the 
exponentially small (at N ^> 1) vicinity of these spacial points. At m> 1 

r c = m 2 , (40) 

and the correlation length diverges when v trends to zero along the special points and there is the Neel ordering in 
this limit. 

At v = 1, when 3\i = — 1, J23 = J13 = 0, the system is divided into the ferromagnetic pairs 1 — 2, 3 — 4, ... . 
Using the second order perturbation theory with respect to [y — 1), we reduce the model (2) at v — ► 1 to the 
effective spin-1 Hamiltonian 

(y — l) 2 M 

Heff = g ^ {~(L n • Ln+i - 1) + L n • L„ +2 

n— 1 

— ' I J n+i)(L„+i • L„ +2 ) - -(L n+ i • L„ +2 )(L„ • L„ + i)| , (41) 

with M spins L = 1. 

The exact singlet ground state wave function of the Hamiltonian (41) can be obtained from the Eqs.(4) and (5). It 
has a form 

*o(M) = P + £ 3 + + - + ^) • • • (^-1 + L+)L+ I - 1, . . . , -1), (42) 

where Lf are raising operators of spin 1. 

The ground state correlation function of (41) K(l,n+ 1) = (Li • L„ + i)/Gm is found from Eqs.(36) and (37). It is 

-n /27m 

Finally, we note that it is possible to calculate higher terms of the perturbation theory and to obtain an effective 
spin-1 Hamiltonians which are proportional to the third, fourth and higher power of the small parameter (y — 1). All 
of them have zero ground state energy as well as (41). 
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FIG. 2. Spin ladder representation of the model (45). Different lines correspond to different exchange interactions. 



III. FRUSTRATED SPIN CHAIN WITH AN ANTIFERROMAGNETIC GROUND STATE. 



A. The model and its exact ground state. 



In the preceding section the spin model at the ferromagnet-antifcrromagnct transition point has been studied. The 
exact singlet ground state wave function at this point is given by Eq.(4). In this section the function (4) will be 
generalized to give the exact ground state wave function of new special one-dimensional frustrated spin-i model. This 
model has the unique singlet ground state (for the cyclic chain) with an exponentially decay of correlations and there 
is a gap to the excitations. 

Let us consider the wave function which depends on two parameters and has the form (4) 

*o(Af) = P *m, (43) 

where 

*M = (st + V X s\ + V 2 s\ ... + V 2 S+)(s+ + V X s\ + V 2 s£ ... + V2Sn) ■ ■ ■ 

x(4-i + ^i4) |1,2,...JV) (44) 

We will construct the Hamiltonian for which ^q(M) is the exact ground state wave function as the sum of the local 
four-sites Hamiltonians 



H — i?l,2,3,4 + #3,4,5,6 + • • • + #JV-3,JV-2,JV-1,JV + #JV-1,JV,1,2 

The Hamiltonians H iti+ i. i+ 2,i+3 are chosen in the form: 



(45) 



-Hi,i+M+2,i+3 — Jl2 

+2Jl 3 



(Sj • Sj+i - -) + (Sj +2 • S i+3 - -) 



(Si • S i+2 - -) + (Si+i • S i+3 - -) 



+2J23(Sj+l • Si + 2 — — ) + 2Ji4(Si • Si + 3 — — ) 



(46) 



Thus, the model has nearest (J12 and J23), next-nearest (J13) and next-next-nearest (J14) neighbor interactions. 
In fact, this model is equivalent to the spin-i ladder with different interactions as it is shown in Fig. 2. 
We demand that ^o(M) 1S the eigenfunction of each local Hamiltonian with the eigenvalue e < 0, i.e. 



ifi,i+i,i+2,i+3*o(M) = e * (M) 
The exchange integrals Jy and the energy e is defined by the Schrodinger equation: 

#1,2,3,4*M = £ 

Let us represent the function $ M in the form 

*M= [qi+v 2 q 2 S + (N -±) + v 2 2 S +2 {N -4)] |1,2,3,4)*m-2, 



(47) 
(48) 
(49) 
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where 



qi = s^sj + v>i(sfs% + s£s£) + insist + "2(1 + ^1)53541 
fa = sf + vis% + (1 + v 2 )s£ + {v\ + 1^2)54 

The equation 

#i,2,3,4?i |1,2,3,4) = £ gi |1, 2,3,4) (50) 

reduces to five equations for four exchange integrals and the energy e. The necessary condition of the existence of a 
solution with e < is 

1 + ^+^2=0 (51) 

So, under this condition there is only one parameter of the Hamiltonian (45). It is convenient to take the value 
/i = 2,X+i as a system parameter. Then, the solution of Eq.(50) at v\ = 2 ^-\ anc ^ ^2 = i^ 2 \ yields the following 
expressions for and e: 

j _ (1-2 M )(2- M ) j _ 1 (2- M ) 2 ( M -l)p 

7, , _ (2m-1)(2-/x)(a»-1) 2 T (2^-l) 2 (^-l) fcol 

J 13 ~ 6(71+1? ' Jl4 - 3(a»+1) 2 ^ > 

_ 3 M (m-1) 2 
£ _ (P+1) 2 



It turns out that the equation 



#1,2,3,492 11,2,3,4) =e q 2 | 1,2, 3, 4) 



with Jij and e given by Eq.(52) is satisfied automatically. 

As it will be proven in Appendix the function with v\ = —v 2 — 1 contains the singlet and the triplet components 
only, i.e. 

S +2 (N)V M = (53) 

Therefore, the last term in Eq.(49) vanishes and \Pm is the eigenfunction of #1,2,3,4 with the eigenvalue e. 

Generally, the Hamiltonian #1,2,3,4 has following eigenstates: one quintet, three triplets and two singlets. Two of 
them (one singlet and one triplet) have the energy e while other four states have higher energies at fi > (fi 7^ 1). 
At /1 < the ground state of #1.2,3.4 is a quintet with zero energy. In Appendix we will show that the wave function 
^ m is the eigenfunction with the eigenvalue e of each local Hamiltonian #, i+ i, i+2 ,i+3 excluding that for i = N — 1. 
Therefore, the ground state energy of the open chain at /i > is (y — l)e and it coincides with the exact lower 
bound of the energy similarly to the model of Section 2. However, in the contrast with the latter the present model 
is four-fold degenerate for the open chain. 

As for the Hamiltonian #jv-i,iv,i,2 the function *m is not its eigenfunction, but 

#at-i,jv,i,2*o(M) = e * (M) 

(see Appendix) and 

N 

#* (M) = —e * (M) 

As it follows from Eq.(52) the spectrum of H(p) coincides with the spectrum of #(/x) which is connected with H(p) 
by transformation 

H(H) = /i 2 #(i) 
M 

This transformation permutes the factors at the third and the last terms in Eq.(46). Therefore, it is sufficient to 
consider the Hamiltonian #(/x) in the region — 1 < fi < 1. 
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The ground state of H is ferromagnetic at fj, < 0. When < ji < 1 the ground state of the cyclic chain is the 
non-degenerate singlet. The point ji = is the ferromagnct-antiferromagnet transition point. At this point the present 
model coincides with the model (2) at the special point v = ^. As it follows from Eq.(52) this transition is the phase 
transition of the first order with respect to ^. 

At /i = | the only non-zero exchange integral is J23 and the ground state consists of non-interacting singlet pairs 
2 — 3,4 — 5, . . . , 1 — N. When /i — > 1 the first term in Eq.(46) dominates and the system is divided into weekly 
interacting ferromagnetic pairs 1 — 2,3 — 4,... . Using the second order perturbation theory with respect to the small 
parameter (fj, — 1) we reduce the Hamiltonian H (Eq.45) to the effective spin-1 model given by 

(M-l) 2 

where 

M 



(54) 



Heff — J]] ^5L n • L„ + i + L„ • L„ + 2 + — (L„ • L„ + i — L n +i • L Il+ 2) 2 — 6 1 , 

and L„ is spin-1 operator. 

The ground state wave function of (54) can be obtained from Eq.(43) at fi = 1. It reads 

vfo(M) = P (1+ - 21 + - . . . - 2L+ )(I + - 2L+ - . . . - 2L+ ) . . . x (L+ _ x - 2L+ )L+ |-1, -1, ... - 1) (55) 

It is remarkable that the function (55) coincides with the ground state wave function of the AKLT model. Therefore, 
the ground state physics of the model given by Eq.(54) and AKLT one is the same, though the Hamiltonians of these 
two models are different. 

B. Spin correlations in the ground state. 

First, we calculate the norm of the ground state wave function. It is convenient to express the function ^> m in 
terms of the parameter /i and to introduce a new function \E"m 



where 



f M = (2M-l)- M *m, 

$ M - [(2m - 1)4 + (2 - n)s+ - (p + 1)S+(N - 2)] [(2 M - l)s+ + (2 - M ) S + - ( M + 1)5+ (A - 4)] 
... [(2 M -1) S +_ 1 + (2- M )s+] |1,2,...A> 

According to Eq. (18) the norm of ^ {M) = P *M is 

G M = (*o(M)* (M)) = \J ®M{y)dy 

For the present model the function <&i{y) is defined by the equation 

$,(!/) = L\y) 1, (56) 



where 

9(1 -m) 2 , (1 + m) 2 



L{y) = 

The solution of Eq.(56) is 



where 

wi =6(m 2 -,u + 1), W2 = 2(/i-2)(2/i-l) 
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This form of reflects the fact that the function 4>m contains the singlet and the triplet components only. 

Thus, Gm is 



Gm = 



Luf 1 + 3wf 



(57) 



As 



< 1 then G M = j^f at M -> oo. 



The spin correlation functions can be found in the similar way as in the Section 2. In analogy with the Eqs.(19)-(20), 
we obtain 



Q(N -2l-l,N) 



-1 j 1 {L M - l -^' l+1 {y))dy, 



(58) 



where 



^l+liv) = (*/+l(PJV-2/-l,iV - l)Po{PN-2l-\,N - 

Carrying out the necessary calculations, we find 

& l+1 (y) = 4(/i + 1) 2 (2 - tfut 1 + 4 +1 + (2/i - l) 2 [(2^ 2 - 2^ + h)^ 1 + (2/i 2 + 10 M - l)^ -1 ] (1 + y) 
Substituting into Eq.(58) we obtain for the spin correlation function K(l, 21 + 2) = K(N — 21 — 1, AT) 



w) = l- »"-"" + 3 ( g ) 

V ' ' 4 2wi 



M-l 



l + 3(t) 



M ' 



#(1,2/ + 2) = - 



3(n 



UJ1LJ2 



If ^-l) 2 (^)' + (^-2) 2 (g) 

/ \ M 



M-l 



where I = 1, 2, . . . M - 1. 

The similar calculations for if(l, 21 + 1) result in 



#(1,2Z + 1) 



3(M+1) ; 



1 + 3(S) 



(59) 



(60) 



M- 



M 



a* 1 + 3 ^y 

The correlators A(2, n) have been obtained by using the symmetry of the system (see Fig. 2). For example, 

K{2,3) = K(1,N), etc. 
In the thermodynamic limit M — > oo and Z <C M Eqs.(59)-(61) reduce to 

3 uj 2 



(61) 



#(1,2) = - 



4wi 



#(l,2Z + 2) 



3(^ + l) 2 (2^-l) 2 /w2 



(62) 



(63) 



#(1,2Z + 1) = fctl)!^ 
2wi V^i, 

The correlators have the exponential decay and the correlation length r c is 

^i(m) 



r c ( M ) = 2 In 1 



W 2 (/i) 



(64) 



(65) 
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It follows from Eq.(65) that the ground state has ultrashort-range correlations. For example, r c (0) = 2 log 1 (|) 
(this value coincides with r c given by Eq.(39) with m = 3) and r c (l) = 21og~ 1 3. In the latter case r c coincides 
with correlation length of the AKLT model. At \i — i the only non-zero correlators are K(2, 3) = if(4, 5) = . . . = 
K(N, 1) = — | in accordance with the dimer character of the ground state. 

The value W2(m) changes the sign at [i = \ and as it follows from Eqs.(62), (63) and (64) the correlators show the 
antifcrromagnetic structure of the ground state at < [i < \ while at \ < yU < 1 there are ferromagnetic correlations 
inside pairs (1, 2), (3, 4), . . . and the antiferromagnetic correlations between the pairs. 



C. Energy gap. 

The Hamiltonian H of the cyclic chain has a singlet-triplet gap A for finite N. It is evident that for ^ = i the gap 
exists for N — > oo and A(i) = i. The existence of the finite gap at the thermodynamic limit in the range < /u, < 1 
follows from the continuity of the function A(^). It is also clear that A(fi) at N — > oo vanishes at the boundary 
points fi — and \i = 1 when the ground state is degenerate and there are low- lying spin wave excitations. 

Unfortunately, the method of the exact calculation of A(/x) in the thermodynamic limit is unknown. For /x — | the 
gap can be found by using the perturbation theory in (fj, — |). In this case A(/i) is 

A( M ) = ± + 0(( M -i) 2 ), M <i 

(66) 

A(M)-i-|(M-|) + 0(( A1 -i) 2 ), M >| 

Eq.(66) shows that A(/z) has a cusp at /j, = |. For the approximate calculation A(/i) we use the trial function of 
the triplet state in the form 

N 

*t = Y^c n 8+* (M), (67) 



where 



n=l 



Ajr 

C2i-i - ae lkl , c 2 i = be ikl , k = —t, t=l,...M 



N 



This trial function gives A(/i) which is 



A ^ = 2 En# m <C m K(n, m) J n , m - 2 Y.n^m \ C l\ K ( n i ™)Jn,m 

I E„ 141 + E„^ m <c m X(n, m) 



N " " ^ ^ 2 2 

is 



Function (68) has minima at k = ^ and fc = 7rfor0</^<i and i < fi < 1, respectively. Then A(/it) at iV oo 

A( A1 ) = §(1 + ^)( M -1) 2 , i<M<l 



3 wi(o)i+w 2 ) ' " 2 

(69) 



The dependance of A(/i) given by Eq.(69) is shown on Fig. 3 together with the results of extrapolations of exact 
finite-chain calculations. Both dependences agree very well for /j, > i. In particular, Eq.(69) correctly reproduces 
(66) at fj, ~ However, A(/x) given by Eq.(69) is not zero at fi = while numerical calculations fit the dependence 
A(/n) ~ y/Jl at — > 0. 

We note that the trial function of the type (67) gives the value 0.7407 for the singlet-triplet gap in the AKLT model. 
This estimate is close to the value 0.7143 obtained by another approach in Ref.20. 

The above consideration refers to the gap in the cyclic chain. The open chain has four-fold degenerate ground state. 
Finite chain calculations show that the spin of the lowest excitation is S — 2 . However, there are also two excited 
singlet and triplet states the energies of which are close to that of S = 2. The difference of these three eigenvalues 
decreases to zero exponentially at — > 00. We expect that these states are degenerate in the thermodynamic limit 
(at /U = I they are degenerate for finite N). The gap in the open chain equals to the difference between the energies 
of the degenerate ground state and the lowest excited one. The extrapolation of the results of finite-chain calculations 
to N — > 00 gives the finite gap in the open chains at < [i < 1. Its value is very close to that for the cyclic chains. 
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FIG. 3. Singlet-triplet gap of the model (45) as a function of the parameter /i. The circles denote the results of the 
extrapolation of an exact finite-chain calculations. The solid line represents the dependence A(/u) given by Eq.(69). 



We have studied a class of the one dimensional quantum spin-i models with competing interactions. The exact 
ground state wave function of these models is found in the special recurrent form. The Hamiltonians of these models 
are the sums of local, non-commuting with each other Hamiltonians; the ground state wave function of the total 
Hamiltonian is the ground state solution for each of them. This means that this ground state wave function is 
superstable 19 with respect to each local Hamiltonian. 

One of the studied models describes the transition point from the ferromagnetic to the spiral state when the energies 
of these two states are equal to each other. It is interesting to compare the exact quantum singlet state with the 
classical one. Both states are the states of a helical type (excluding some special cases) and their spin correlation 
functions are identical in the thermodynamic limit though quantum effects are essential for finite chains. This fact is 
rather surprising for the one-dimensional model with spin |. 

Another exactly solvable Hamiltonian has the antifcrromagnctic ground state. This state is non-degenerate for the 
closed chains and is four-fold degenerate for the open ones. The Hamiltonian depends on the one parameter /j, and 
there are two special values, /j, — and /i = 1 where the singlet and the ferromagnetic states are degenerate. The 
value fi — is the ferromagnet-antiferromagnet transition point where the phase transition of the first order with 
respect to [i occurs. 

The ground state is characterized by the exponential decay of correlators with a very short correlation length, and 
there is the gap in the excitation spectrum at < /j, < 1. Thus, this model has all properties suggested by Haldane 16 
for the one-dimensional Heiscnbcrg antifcrromagnet with integer spin. The first model for which all these properties 
have been proved rigorously is the AKLT model. Our model is the one with spin i. Affleck and Lieb 21 have shown 
for the translationally invariant and the isotropic Heisenberg Hamiltonians that for half-integer spin chain either the 
excitation spectrum is gapless or the ground state is degenerate. The existence of the finite gap in our model does not 
contradict to Afflcck-Lieb theorem because this model is not translationally invariant. It has two sites in the unit cell 
and is equivalent to the special ladder model. Moreover, in the limit ^ — > 1 its ground state wave function reduces to 
that for the AKLT model. 
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VI. APPENDIX 

We prove that square of the total raising operator annihilates ^ m if the condition (51) is satisfied. The recurrent 
equation for S +2 (N)^ M is 

S+ 2 (N)<f M = 2(1 + v x + u 2 )stst\l, 2) S+(N - 2)* M -i+ 

(Al) 

+ [(1 + 2v 2 )s+ + fa + 2u 2 )s+ + u 2 S+{N - 2)] |1, 2) S+ 2 {N - 2)* M -i 
The first term in Eq.(Al) vanishes under the condition (51) and, therefore, 
S +2 (N)V M = [(1 + 2u 2 )s+ + fa + 2v 2 )s+ + v 2 S+(N - 2)] 

x [(1 + 2v 2 )s+ + fa + 2v 2 )s+ + v 2 S+{N - 4)] ... x (4^ + s+ f [s+_ x + u lS +] |1, 2, . . . N) = 

This equation means that the wave function contains the singlet and triplet components only. 
Now we prove that &o(M) is the eigenfunction of each local Hamiltonian in (45). Of course, consequentely the 
same will be true for 5 , (M). 

The function \P m satisfies the recurrent equation 

*m = (u 12 $ M -i-(m+1)S + (^-2)$m-i) |1,2), (A2) 

where 

u 12 = {2fx - l)s+ + (2 - ^)s+ 

Let us consider functions ip^ = S + (N)^ M , V M = S~(N)y M and \m = S~ (N)S + (TV)* M • The recurrent 
equations for these functions are obtained from (A2) using (Al). They are 

Vm = ((m + 1)44*m-i+£i2¥>m_i) |1,2), 
<Pm = ((m + 1)(*m-i-Xm-i)+ Wm-i) I 1 ' 2 )' ( A3 ) 

XM= ("12 - Wl 2 )*M-l + V12XM-I + (/" + l)(s+ s 2>M-l - <4f-l) |1,2), 

where 

«i 2 = -(2-/i)4-(2/i-l)4 
Eqs.(A2)-(A3) can be written in a matrix form 

R(M) = D 12 R(M - 1), (A4) 

where R(M) and D 12 are (2 x 2) matrices 

R(M) - f * M ~ XM ^ 



Dl2 -{-{?+ 1) -12 J !l ' 2) 



Therefore, i?(M) is 

R(M) = D 12 D 34 . . . D N _ hN (A5) 
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As ^ m contains the singlet and triplet components only, the projection of ^ m onto the singlet is 

*o(M) - P *m = 2* M - Xm (A6) 

It follows from Eqs.(A5)-(A6) that 

*o(M) - Tr(D 12 D 34 . . . D N _ hN ) (A7) 

This form of 'I'o(Af) is similar to the matrix product wave function of the AKLT model and its generalizations 
which has been found in Ref.15. Each of four matrix elements of R(M) is the eigenfunction of the local Hamiltonian 
#i,i+i,i+2,i+3 for i = 1, 3, ... N — 3 because of matrix elements of the product D iti+ iD i+ i.i + 2 are the eigenfunctions 
of this Hamiltonian. Besides, it can be proved 15 that the four matrix elements of R(M) are the only ground states of 
(45) and, therefore, the ground state of the open chain is four-fold degenerate. 

It is easily to check 15 that the triplet wave functions R\i{M) and R2i(M) are not eigenfunctions of fljv-i,jv,i,2- 
On the other hand, using cyclic permutations of matrices under the trace, we have 

iJjV-l,JV,l,2*o(M) = e tf (M) 

and, therefore, the ground state of the cyclic chain is the non-degenerate singlet. 
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